Mon. Not. R. Astron. Soc. 000,[T][8](2003) Printed 2 February 2008 (MN MbX style file v2.2) 



Variable stars in the Magellanic Clouds: 
II. The data and infrared properties 

Yoshifusa Ita^'^*t, Toshihiko Tanabe^ Noriyuki Matsunaga^ Yasushi Nakajima^, Chie 
Nagashima^, Takahiro Nagayama"^, Daisuke Kato"^, Mikio Kurita"^, Tetsuya Nagata"^, 
Shuji Sato"^, Motohide Tamura^, Hidehiko Nakaya^ and Yoshikazu Nakada^'^ 

^Institute of Astronomy, School of Science, The University of Tokyo, Mitaka, Tokyo 181-0015, Japan 

'^Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency, Sagamihara, Kanagawa 229-8510, Japan 
^National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, Japan 
^Department of Astrophysics, Nagoya University, Chikusa-ku, Nagoya 464-8602, Japan 

^Subaru Telescope, National Astronomical Obsen'atory of Japan, 650 North A'ohoku Place, Hilo, HI 96720, U.S.A. 
^Kiso Observatory, School of Science, The University of Tokyo, Mitake, Kiso, Nagano 397-0101, Japan 



Received - / Accepted - 



ABSTRACT 

The data of 8,852 and 2,927 variable stars detected by the OGLE survey in the Large and 
Small Magellanic Clouds are presented. They are cross-identified with the SIRIUS JHK sur- 
vey data, and their infrared properties are discussed. Variable red giants are well separated on 
the period-7 — K plane, suggesting that it could be a good tool to distinguish their pulsation 
mode and type. 

Key words: galaxies: Magellanic Clouds - stellar content - stars: AGB and post-AGB - 
variables - infrared: stars - surveys 



■ 1 INTRODUCTION 

■ Because of their proximity and comparatively well-known dis- 
tances, the Large and Small Magellanic Clouds (LMC and SMC, 
respectively) are the ideal places to study stellar evolution. To 
search for the gravitational lensing events, several projects are now 
underway monitoring the sky toward the Magellanic Clouds and 
also the Galactic bulge, where there are a large number of back- 
ground stars that are potential targets for microlensing. As the nat- 
ural by-products of these surveys, many variable stars were newly 
found. The Optical Gravitational Lensing Experiment (OGLE) 
project is one of such s urveys and it s / band time-series data in 
the Magellanic Clouds (Zebruhetal. 2001) obtained during four 
years of OGLE-II run lUdalski, Kubiak & Szvmanski 1997; Jan- 
uary 1997 - November 2000) is now available over the Internet 
(OGLE homepage; http://sirius.astrouw.edu.pl/'^ogle/ 1. 

To detect dust-enshrouded variables that have possibly es- 
caped from the ongoing/previous optical projects, we started the 
JHK monitoring survey toward the Large and Small Magellanic 
Clouds using the InfraRed Survey Facility (IRSF) at the South 
African Astronomical Observatory Sutherland station. The IRSF 
consists of a dedicated 1.4m telescope to which is attached a near 
infrared (NIR) camera SIRIUS, which can observe the sky in the 
three wave bands (JHKs) simultaneously. The filter characteristics 
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Table 1. The filter characteristics. 





J 


H 




Central wavelength [um] 


1.25 


1.63 


2.14 


Bandwidth |jum] 


0.17 


0.30 


0.32 



are shown in table Q The field of view of this system is about 
1.1'xl.l' with a scale of 0.453"/pixel. Details of the instrument 
are found in Nagashima et al. ( 1999) and Naaavarna et al. (200^. 
This campaign was initiated in December 2000, and so far, many 
dust-enshrouded variables were found and their light variation data 
has been steadily accumulating (see Ita et al. 2002 for the first re- 
sults). 

In the previous paper llta et all200'4l hereafter referred to as 
Paper I), we cross-identified the OGLE data and the single-epoch 
SIRIUS JHK data from the monitoring survey, and studied the pul- 
sation properties and metallicity effects on period-/^ magnitude re- 
lations by comparing the variable stars in the LMC and SMC. In 
this paper, we present the data of the variable stars in the Magel- 
lanic Clouds and discuss their infrared properties. 



2 DATA 

The total areas covered by the SIRIUS variable star survey are 
three and one square degrees in the LMC and SMC, respectively. 
These areas are divided into 27(9) regions of 20'x20' meshes in 
the LMC(SMC). Tables |2l and |3] define the central coordinates of 
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Table 2. Survey regions in the Large Magellanic Cloud 



Name Right Ascension Declination 



Equinox: J200() 



LMC0507-6840 


OS*" 




00.28" 


-68° 


40' 


00.00" 


LMC0507-6900 


OS*" 


06"" 


50.29" 


-69° 


00' 


00.00" 


LMC0507-6920 


OS*" 


06" 


39.97" 


-69° 


20' 


00.00" 


LMC05 10-6920 


OS*" 


10"' 


26.65" 


-69° 


20' 


00.00" 


LMC05 10-6940 


OS*" 


10"' 


19.53" 


-69° 


40' 


00.00" 


LMC05 11-6900 


OS*" 


10"' 


33.53" 


-69° 


00' 


00.00" 


LMC05 14-6900 


OS*" 


14m 


16.75" 


-69° 


00' 


00.00" 


LMC05 14-6920 


OS*" 


14m 


13.31" 


-69° 


20' 


00.00" 


LMC05 14-6940 


OS*" 


14"' 


09.77" 


-69° 


40' 


00.00" 


LMC05 18-6900 


OS*" 


18"' 


00.00" 


-69° 


00' 


00.00" 


LMC05 18-6920 


OS*" 


18"' 


00.00" 


-69° 


20' 


00.00" 


LMC05 18-6940 


OS*" 


18"' 


00.00" 


-69° 


40' 


00.00" 


LMC0522-6920 


OS*" 


21"' 


46.67" 


-69° 


20' 


00.00" 


LMC0522-6940 


OS*" 


21"' 


50.23" 


-69° 


40' 


00.00" 


LMC0522-7000 


OS*" 


21"' 


53.89" 


-70° 


00' 


00.00" 


LMC0526-6920 


OS*" 


25™ 


33.35" 


-69° 


20' 


00.00" 


LMC0526-6940 


OS*" 


25™ 


40.45" 


-69° 


40' 


00.00" 


LMC0526-7000 


OS*" 


25™ 


47.81" 


-70° 


00' 


00.00" 


LMC0529-6920 


OS*" 


29™ 


20.01" 


-69° 


20' 


00.00" 


LMC0530-6940 


OS*" 


29™ 


30.67" 


-69° 


40' 


00.00" 


LMC0530-7000 


05" 


29™ 


41.71" 


-70° 


00' 


00.00" 


LMC0533-6940 


05" 


33™ 


20.91" 


-69° 


40' 


00.00" 


LMC0534-7000 


05" 


33™ 


35.61" 


-70° 


00' 


00.00" 


LMC0534-7020 


05" 


33™ 


50.83" 


-70° 


20' 


00.00" 


LMC0537-6940 


05" 


37™ 


11.13" 


-69° 


40' 


00.00" 


LMC0537-7000 


05" 


37™ 


29.52" 


-70° 


00' 


00.00" 


LMC0538-7020 


05" 


37™ 


48.53" 


-70° 


20' 


00.00" 



Table 3. Survey regions in the Small Magellanic Cloud 



Name Right Ascension Declination 









Equinox 


: J2000 




SMC0050-7250 


00" 


50™ 


' 28.95" 


-72° 


50' 


00.00" 


SMC0050-7310 


00" 


50™ 


' 23.74" 


-73° 


10' 


00.00" 


SMC005 1-7230 


00" 


50™ 


' 33.95" 


-72° 


30' 


00.00" 


SMC0055-7230 


00" 


55™ 


' 00.00" 


-72° 


30' 


00.00" 


SMC0055-7250 


00" 


55™ 


' 00.00" 


-72° 


50' 


00.00" 


SMC0055-7310 


00" 


55™ 


' 00.00" 


-73° 


10' 


00.00" 


SMC0059-7230 


00" 


59™ 


' 26.04" 


-72° 


30' 


00.00" 


SMCOlOO-7250 


00" 


59™ 


' 31.04" 


-72° 


50' 


00.00" 


SMCOlOO-7310 


00" 


59™ 


' 36.25" 


-73° 


10' 


00.00" 



them. Each of the 20'x20' regions is further subdivided into nine 
sub regions due to the field of view of the SIRIUS, which is about 
7.7'x7.7'. 

We cross-identified OGLE variables in the Magellanic Clouds 
jZebruh et al. 2001) with the SIRIUS NIR sources that have been 
detected at least two of the three (JHK) wave bands, and then 
determined their pulsation perio ds by using the Ph ase Dispersion 
Minimization (PDM) technique <StelIingweril978h . Details of the 
cross-identifications and period findings are described in Paper I. 
We did not analyze multi-periodic stars (e.g., Beddine et al. 1998) 
and/or variable stars that show too complex light curves in the 
phase space (see the bottom panel of figures|2|and |3j to find a pre- 
dominant period. We made two catalogs presenting data of variable 
stars in the LMC and SMC. They contain the following informa- 
tion: OGLE name, pulsation period in days, statistical parameter 6 
calculated by the PDM, pulsation amplitude (Al = /^ax — ^min)> in- 
tensity mean / magnitude < / >, JHK magnitudes and positional 
difference r between the SIRIUS coordinate and the OGLE coor- 
dinate in arcsec. The is a measure of the regularity of the light 



variation, being near for the regular variation and near 1 for the 
irregular variation. 

Table |4| is a sample that shows the first five records of the 
LMC catalog. The full version of the data, including the SMC cat- 
alog is available in the on-line version of this paper. In cases the 
JHK measurements are saturated or below the detection limit, the 
corresponding columns contain the value 99.999. The JHK mag- 
nitudes in the catalog are referred to the Las Campanas Observa- 
tory (LCO) system ones b ased on observations of a few dozen stars 
from lPersson et all il998h and are not corrected for the interstellar 
absorption. We used the following colour equations 

/ -0.007 ±0.002 \ 
^LCO=^+ -0.005±0.002 + (1) 
\ -0.003 ±0.002 / 

0.016±0.011(y-//) \ 
0.009 ±0.005(y-/f) 
0.000 ±0.000(tf-/f) / 

/ 0.001 ±0.002 \ 
Hlco=H+1 -0.002 ±0.002 + (2) 
\ 0.000 ±0.002 / 

-0.001 ±0.000(/-//) \ 

0.004 ±o.ooo(y-/f) 

0.004±0.000(/f-A:) / 

/ -0.001 ±0.004 \ 
Klco = K+\ -0.001 ±0.004 + (3) 
\ 0.005 ±0.004 / 

-0.002 ±0.000(/-//) \ 

-0.002 ±o.ooo(y-/f) 

-0.041 ±0.004(/f-/f) / 

to transform the IRSF instrumental magnitudes to the LCO system 
based ones. The J — K ones were usually used, but J — H(H — K) 
ones were employed if the K( J) magnitudes were unavailable. Note 
that the SIRIUS measures the K^ magnitudes and they were con- 
verted to the LCO K magnitudes. 

To enhance the image quality, we stacked the top 10 best see- 
ing images' (typically the seeing ranges from about 0.9" to 1.2" 
at Ks) that we have been collecting in the ongoi ng SIRIUS mon- 
itoring survey. We used the ISIS.V2.1 package iAIard & Luptonl 
Q998; Alard 2000) to build a composite image that comprises 100 
dithered 5 sec exposures. Photometry was performed on the com- 
posite image with DoPHOT i Schechter. Mateo & SahJl993h . and 
then, the measured magnitudes were calibrated with the single- 
epoch photometry data by matching the magnitudes of the non- 
variable stars. Eventually, the photometry has a signal-to-noise ra- 
tio (S/N) of about 10 at 17.40, 17.38 and 16.63 mag in J, H and K, 
respectively. Also, traditional experiments of adding artificial stars 
to the composite images revealed that the 90% completeness (de- 
tection probability) limits at J, H and K are about 17.25, 17.22 and 
16.45 mag, respectively Llta.2004bi) . 

3 DISCUSSION 

In the following discussions, we use the JHK magnitudes that are 
corrected for the interstellar absorption based on the relations in 

' These 10 images were taken on different days and hence under the differ- 
ent conditions, and a single image comprises 10 dithered 5 sec exposures. 
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Table 4. The first 5 records in the LMC catalog. This is a sample of the full version, which is available in the electronic version of Monthly Notices. 



Name 


Period 


e 


AI 


</> 


J 


H 


K 


r 




[days] 








[mag] 






[ai'csec] 


OGLE050437.94-692217.2 


0.460 


0.178 


1.025 


19.250 


99.999 


16.056 


16.030 


2.260 


OGLE050438.43-691527.2 


49.892 


0.837 


0.089 


13.951 


13.559 


10.786 


11.337 


0.306 


OGLE050440.26-693039.7 


42.783 


0.866 


0.160 


14.189 


12.786 


11.880 


11.599 


0.466 


OGLE050440.29-692708.3 


39.200 


0.661 


0.096 


14.734 


13.423 


n.6\l 


12.333 


0.465 


OGLE050440.76-692416.0 


96.365 


0.358 


0.317 


14.060 


12.294 


11.491 


11.063 


0.574 




Table 5. Colour mapping table 



Log P [days] 

Figure 1. Variable stars in the LMC (upper panel) and SMC (lower panel) 
are classified into several groups according to their location on the period-A" 
magnitude plane. 



lKoornneed<1982h . assuming R = 3.2. We adopted (Ay, Ah, Ak) ~ 
(0.129, 0.084, 0.040) and (0.082, 0.053, 0.025), coiresponding to 
the total mean reddening of Eg^y = 0.137 and 0.087 as derived by 
[Udalski et al. ( 1999) for the OGLE's observing fields in the LMC 
and SMC, respectively. The reddeni ng of the Magellanic Clouds 
is still somewhat controversial (e.g.. IWesterlund . 1997) . and it is 
certain that this type of statistical average approach can not account 
for individual corrections. However, the adopted reddening should 
not affect any of the conclusions of this paper. The coefficients we 
calculated in equations|4|toQdepend on reddening, but the effects 
should be small. 



3.1 Classification of the variable stars on the period-A" 
magnitude plane 

IWood et'all {l99'f) discovered that variable red giants in the LMC 
form parallel sequences in the period-A" magnitude plane and sug- 
gested that these seq uences can b e interpreted as the differences in 
the pulsation modes jWooJ2000l) . 

In Paper I, we classified variable stars into nine (LMC) and 
eight (SMC) prominent groups based on their locations on the 



Label 


Colour 


Population 


A- 


Cyan 


RGB variables & metal poor and old AGB variables 


A+ 


Yellow 


less regularly pulsating AGB variables 


B- 


Orange 


RGB variables & metal poor and old AGB variables 


B+ 


Green 


less regularly pulsating AGB variables 


C 


Blue 


Mira variables pulsating in the 1 st overtone mode 


C 


Red 


Mira variables pulsating in the fundamental mode 


D 


Steel-blue 


Some obscured variables & Unknown variables 


F 


Magenta 


Cepheid variables pulsating in the fundamental mode 


G 


Purple 


Cepheid variables pulsating in the 1st overtone mode 



period-A^ magnitude plane. Refer to Paper I for the details of the 
grouping. For the easier reading and completeness of this paper, 
we show the period-AT magnitude diagram with classification boxes 
in figureQ and also for conciseness, we tabulated the colour map- 
ping chart in table|5| Stars outside the classification boxes are rep- 
resented by black tiny dots. Throughout this paper, we use the same 
labels and colours for each group. In the LMC, there is a clump of 
variables around logP < —0.1 and K > 15.0. This is a group of RR 
Lyrae variables, which we do not analyze in this paper. 

Roughly speaking, five types of variable stars constitute the 
majority of the catalog; Cepheids, Miras, Semi-regulars, Irregular 
variables and Eclipsing binaries. However, telling Semi-regulars 
from Mira variables is rather difficult and there are not necessar- 
ily dee p physical criteria. Several authors (e.g., Alard et al."20()]i 
ICioni et al. 2001J use a simple criterion that is based on a certain 
limiting pul sation amplitude t o sepa rate Semi-regulars and Miras. 
Meanwhile, iKiss & Beddind j2003l 2004) showed that there is a 
good correlation between the amplitude and mode of pulsation. In 
this paper, we do not separate Semi-regulars and Miras explicitly. 
Instead, we distinguish "regularly pulsating variables, (6 ^ 0.55)" 
from "less regularly pulsating variables, (6 > 0.55)" relying on the 
parameter "0". In figures|2|and|3| we show the representative light 
curves of two variable stars on sequence C, having similar periods 
(~109 and ~126 days) and the same amplitudes (A/ ~0.490 mag) 
but different 6s. 



3.2 Colour-magnitude and colour-colour diagram 

The NIR data from the SIRIUS survey enabled us to study the 
infrared properties of the OGLE variables. Figure |4] and |5| show 
the NIR colour-magnitude and colour-colour diagram of the OGLE 
variables. In the inset of the figure|4] we show the colour-magnitude 
diagram of the SIRIUS data (including non-variables), closing-up 
around the tip of the first giant branch (TRGB). By comparing the 
insets and the main diagrams, it is clear that the stars of group 
and accumulate at the exact l ocation of the T RGB l lta et all 
2002; Kiss & Beddina 20o| 2004: llta et all2004j) . 

Most of the stars from group D are found at the place a 
bit brighter than the TRGB. Wood et al. 1 1999j i dentifi ed the se- 
quence D with suspected binaries, and lWood et aljj2004h discussed 
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Figure 2. The light curve, theta spectrum and phased hght curve (from top 
to the bottom) of the star OGLE051707.45-692855.3. 
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Figure 3. The same as frgurelU but for the star OGLE052353. 10-694718.4 
to show a similar star with different value of 6. 



the possible explanations for the D variables. On the other hand, 
some of the D variables are very red (J — H > 1.8, H — K > 
1.5 and y-A: > 3.2) and relatively bright (11.0 < A: < 12.5), 

i iroperties associated with c arbon-rich stars and dusty AGB stars 
Nikolaev & Weinber J2o'ool) . In the present sample, 14 out of 471 
(6 out of 71) D stars in the LMC (SMC) have these properties. A 
possible explanation for these red stars in the group D is that, these 
stars could in fact be very dusty AGB stars that ar e subjected to se- 
vere c ircumstellar extinctions, even in the K band. lWhitelock et alJ 
i2003t ) studied obscured variables in the LMC and showed that the 
K magnitudes of obscured stars cover a wide range at a given pul- 
sation period. For instance, their data indicates that the difference 
can be more than two magnitudes in the case of the carbon-rich 
stars. It is interesting to note that their red stars (those with long 
periods and K magnitudes fainter than ~11 mag) clearly locate the 
brighter part of se quence D on the period-A^ magnitude plane. In 
addition. |Kiss & Beddinj l2003l) suggested that there is a group of 
stars located below the Mira sequence (sequence C) and identified 
them with dusty AGB stars. 



3.3 Period-colour diagram 

Figure|6|shows the J — K colours of OGLE variables in the Magel- 
lanic clouds as a function of their periods. The insets are the close- 
ups for the sake of clarity. A general trend that the colours get red- 
der with increasing period can be seen. 



The dashed lines are the least-square fit of a linear relation 
to the Cepheid variables (F and G stars), which yields the period- 
J — K colour relations of 

y-/i:[mag] =0.147(±0.010)log/'[days]+0.248(±0.006) (4) 

for the LMC (o = 0.062) and 

y-/i:[mag] =0.167(±0.009)logP[days]+0.272(±0.004) (5) 

for the SMC (<J = 0.074), respectively. The quoted errors are the 
la errors of each coefficient. Although the fundamental (group F) 
and first overtone (group G) Cepheids are clearly separated in the 
period-A^ magnitude plane, we don't see any significant difference 
between them in this diagram. This means, if we compare the fun- 
damental and first overtone Cepheids with the same pulsation peri- 
ods, the stellar radii of first overtone Cepheids should be larger than 
those of fundamental Cepheids because their stellar temperatures 
are nearly the same but the first overtone Cepheids are brighter than 
the funda mental one s. Theories predict exactly the same conclu- 
sion (e.g. jBono et all2 001). and if the stellar radius of the Cepheid 
variables could be measured accurately, it will help to identify the 
pulsation modes. 

The less regularly pulsating red giants (group and B^) also 
follow a tight sequence on the period-/ — K colour plane. The thick 
solid lines in the figure indicate the least-square fit of a linear re- 
lation to the stars from the four groups, whose J — K colours and 
periods are within the ranges, 0.7 < J — K < 1.3 and 1.2 < logP < 
2.1. Their period-y — K colour relations are also calculated and they 




J - K [mag] 



Figure 4. NIR colour-magnitude diagram of OGLE variables in the LMC 
(upper panel) and the SMC (lower panel). The inset shows the plot of the 
all SIRIUS data (including non-variables) around the A^-band luminosity of 
the TRGB. 

are, 

y-/i:[mag] =0.231 (±0.004) log P[days]+ 0.671 (±0.007). (6) 
for the LMC (o = 0.058) and 
y-/i:[mag] =0.272(±0.011)log/'[days]+0.560(±0.019). (7) 

for the SMC (o = 0.068), respectively. 

It is obvious that most of the C and C' stars in the LMC, and 
also some S+ stars and most of C stars in the SMC do not follow the 
extension of these relations. In Paper I, we confirmed the conclu- 
sion o f Wood & Seb o ( 1996) and suggested that stars on sequences 
C and C' are Mira variables pulsating in the fundamental and first 
overtone mode respectively, by comparing their model calculations 
with our observational data. Figure |6| clearly shows that C and C' 
stars behave very similarly on the period-y — K colour plane. This 
corroborates the idea that the C and C' stars are the same type of 
stars (i.e., Miras), but only different in the pulsation modes. Also, it 
is likely that the "some 5+ stars" in the SMC are the counterparts 
of the C' stars in the LMC. 

In earlier work . Feast et al. ( 1989), iGlass et al] <1995h and 
IWhitelock. Mararig" &: Feast ( .2000.) obtained period-7 — K colour 
relations of Mira variables in the LMC, Sgr I Baade window of the 
Galactic Bulge and solar neighborhood, respectively. These rela- 
tions are all based on the S AAO system. The atmospheres of Mira 
variables are rather complex, making their quantitative compar- 
isons difficult between different filter systems. Thus here we make 
only qualitative discussions. Just to get a rough idea, we indicated 
in the upper panel of figure|6|the period-y — K colour relations ob- 
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- 
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H - K [mag] 

Figure 5. NIR colour-colour diagram of OGLE variables in the LMC (upper 
panel) and the SMC (lower panel). 



tained by iFeast et all i 19891) for carbon- (labe led I) and oxygen- 
rich (labeled II) Miras and bv lWhitelock et alj i2000) (labeled III) 
after referring to the LCO system by assuming (7 — A^)lco = 
0.928(y-A:)sAAO -0.004 I Carpenter 2001). Of course we can not 
tell exactly which of our stars are carbon- or oxygen-rich, but it 
is statistically fair to say that stars with J ~ K co lours redder than 
1.4 are primarily carbon-rich stars iNikolaev &_We inberg 2000). In 
this point of view, we confirmed the result of lFeast et alJil989 ) that 
oxygen- and carbon-rich Miras follow different period-7 — K colour 
relations. 

Now that we see the C and C' stars separate on the period- 
J — K colour plane, it will be interesting to search for the Galac- 
tic counterpart of the first overtone Miras (C' stars) by using this 
tool. Then we can redetermine their distances using their exclu- 
sive period-A^ magnitude relation, which we obtained in Paper I. 
It will help to understand the differences between the fundamen- 
tal and first overtone Miras through more detailed studies that are 
technically difficult in the Magellanic Clouds (i.e., spectroscopic 
work and mid-infrared observations etc.). Because not only C and 
C' stars but also the other variables separate on the period-y — K 
colour plane, the same argument is true for them, such as D stars. 

However, before applying this tool to different environments, 
we should know whether it depends on chemical abundance. 
IWrav. Ever & PaczvnskJ i2004h showed that small amplitude red 
giant variables in the Galactic Bar can be separated into two groups, 
and they follow different period-colour relations. Also, they suggest 
that their A and B stars correspond to type and variables, re- 
spectively of Paper I. The insets of the figure|6|clearly confirm their 
result, showing that A^ and B^ are well separated on the period- 
J — K color plane. This infers that the period-y — K relations we 
derived in equations |6| and are merely the general trends of the 
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Figure 6. Period-colour diagram of OGLE variables in the LMC (upper 
panel) and the SMC (lower panel). The dashed and thick solid lines in- 
dicate the least-square fits of linear relation to each distinct sequence. 
The thin solid lines labeled by I to III in the uppe r panel are the period- 
J — K colour relations f or Mira variables obtained bv lFeast et alj ^989) and 
IWhitelock et alJi200Ch (see text). The insets show the close-ups of less reg- 
ularly pulsating variables. 



variable red giants, and in fact there are exclusive period-colour re- 
lations for each population. Reikubaetal. 1 2003) studied infrared 
properties of long period variables in Cen A. Although they did 
not mention it, it is likely that there are two groups in their period- 
colour diagrams (see their figure 10). We suggest that their short- 
period (logP ~ 2.45) and long-period (logP ~ 2.65) groups could 
correspond to type C' and C stars in the LMC, respectively. These 
facts strengthen the idea that, period-colour diagram is a good and 
universal tool, at least to tell the type of variable stars within a sys- 
tem. 

Leb zelter. Schultheis & MelchioJ <2002h studied the AGAP- 
EROS/DENIS variables in the LMC and showed that so-called 
Semi-regular variables follow a period-colour relation but the "reg- 
ular" variables do not (see their figure 12). To see the regularity 
of the light variation quantitatively, we show the relationship be- 
tween the statistic parameter 6 and period in figure Q The figure 
clearly confirmed the Lebzelter's suggestion and shows that stars 
that follow period-y — K relation (A* and stars with logP < 
2. 1) are less regularly pulsating, and those that don't (S+ stars with 
logP > 2. 1, C' and C stars) are regularly pulsating. 



3.4 Period-amplitude (A/) and colour-amplitude diagram 

Figure |8| is a plot to show the relationship between the pulsation 
period and amplitude (A/). The pulsation amplitude is determined 
from the OGLE's / band data, given as A/ = /^ax ^^min- Note that 
all of the original OGLE light curves were carefully eye-inspected. 



SB 




Figure 7. Period-9 diagram of OGLE variables in the LMC (upper panel) 
and the SMC (lower panel). 



and we calculated the A/ after eliminating the obvious spurious data 
points. 

The figure shows that the fundamental Cepheids show larger 
amplitude variations than the first overtone Cepheids. Also, the up- 
per envelope of Cepheid amplitude goes down with increasing pe- 
riod. This trend contrasts markedly with that of red giants, because 
their pulsation amplitudes generally go up with increasing period. 

Variables stars with K magnitudes below the TRGB (A^ and 
stars) pulsate with very small amplitude. With a few excep- 
tions, most of the stars from D group also pulsate with relatively 
small amplitude. The large amplitude pulsation (A/ > 0.9) seems 
to occur only among the long-period (logP > 2.3) variables. If one 
had to separate Semi-regular from Mira variables based on pulsa- 
tion amplitudes, the threshold would be about A/ ~ 0.9, because 
one can see a discernible gap around there. 

Figure|9|shows the relationship between the J — K colour and 
pulsation amplitude (A/). The variable red giants show the clear 
trend t hat the amplitude gets l a rger w ith increasing colour. Accord- 
ing to iNikolaev & Weinberg! yOOu), obscured carbon-rich AGB 
stars have the J — K colours redder than about 2.0 mag. Almost 
all of such dusty carbon-rich stars in the Magellanic Clouds have 
large pulsation amplitudes. Some of the relatively blue stars (0.8 
< 7 — A" < 1.2) on the sequence C also show large amplitude pul- 
sation. Judging from their J — K colour and being on the sequence 
C, they are likely to be oxygen-rich Mira variables. It is remarkable 
that the carbon-rich Miras tend to have greater /-band amplitudes 
at redder J — K colour, but the amplitudes of oxygen-rich Miras are 
almost independent of the colour. This might reflect the difference 
in the molecular opacities at /-band between the two groups. 
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Figure 8. Period-amplitude(A/) diagram of OGLE variables in tlie LMC 
(upper panel) and the SMC (lower panel). 




Figure 9. Colour-amplitude(A/) diagram of OGLE variables in the LMC 
(upper panel) and the SMC (lower panel). 



3.5 Similarities and differences between the LMC and SMC 
samples 

In a large sense, variable stars in the Magellanic Clouds behave 
similarly in the figures we discussed so far. However, on closer 
view, there are differences between variable stars in the two galax- 
ies. In the colour-amplitude diagram (figure |9j, we see smaller 
number of oxygen-rich large amplitude Mira variables in the SMC 
(counterpart of LMC Miras with J — K colours bluer than 1 .4) com- 
pared to the LMC. This is probably because the difference in aver- 
age chemical abundance and/or in age between the tw o galaxies 
(e.g.. lCioni & Habin j20o'3LlMouhcine & Lanconl20ol . 



4 SUMMARY 

We presented the data of 8,852 and 2,927 variable stars in the Large 
and Small Magellanic Clouds. Based on these data, we discussed 
the infrared properties of the variable stars. We showed that period- 
J — K colour diagram is a good tool to tell the pulsation modes of 
the Mira variables. A follow-up campaign searching for the Galac- 
tic counterpart of the first overtone mode Miras would tell us the 
difference between them and fundamental mode Miras. 
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